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ABSTRACT 

Wide-separation lensed QSOs measure the mass function and evolution of massive 
galaxy clusters, in a similar way to the cluster mass function deduced from X-ray- 
selected samples or statistical measurements of the Sunyaev-Zeldovich effect. We com- 
pute probabilities of strong lensing of QSOs by galaxy clusters in dark energy cosmolo- 
gies using semianalytical modelling and explore the sensitivity of the method to various 
input parameters and assumptions. We highlight the importance of considering both 
the variation of halo properties with mass, redshift and cosmology and the effect of 
cosmic scatter in halo concentration. We then investigate the extent to which observa- 
tional surveys for wide-separation lensed QSOs may be used to measure cosmological 
parameters such as the fractional matter density 17m, the rms linear density fluctu- 
ation in spheres of Sh^Mpc, <7g, and the dark energy equation of state parameter 
w. 

We find that wide-separation lensed QSOs can measure as and in an equiv- 
alent manner to other methods such as cluster abundance studies and cosmic shear 
measurements. In assessing whether lensing statistics can distinguish between values of 
w, we conclude that at present the uncertainty in the calibration of as in quintessence 
models dominates the conclusions reached. Nonetheless, lensing searches based on cur- 
rent QSO surveys such as the Two-degree Field and the Sloan Digital Sky Survey with 
10 4 -10 5 QSOs should detect systems with angular separations A8 > 5" and hence can 
provide an important test of the standard cosmological model that is complementary 
to measurements of cosmic microwave background anisotropies. 
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1 INTRODUCTION 

The first year Wilkins on Microwave Anis otropy Probe 
(WMAP) measurements (ISpergel et al.ll2u03ft of the cosmic 
microwave background anisotropy seem to confirm that the 
universe is flat and that there is a sufficient density of "dark 
energy" to accelerate the expansion of the universe. It is not 
yet clear whether the dark energy is a cosmological constant 
which is spatially uniform and constant with time, or in- 
stead a dynamic cosmic fiefd, so-called quintessence, which 
varies in space and time. A combination of the WMAP ob- 
servations with other astronomical datasets has provided a 
good measurement of the dark energy equati on of state pa- 
rame ter w < —0.78 (95% confidence limit) dSpergel et alJ 
2003). However, the combination of datasets in this way 
means that independent tests of both the values of cosmo- 
logical parameters and the underlying assumptions cannot 
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be made. Statistical tests based on the abundance of massive 
structures can be an alternative test of the standard cosmo- 
logical model in which structures have grown with the ex- 
pansion of the Universe from initially Gaussian primordial 
fluctuations. Observational signatures of clusters at z ^ 3 
contrast with observations of the cosmic microwave back- 
ground anisotropies at z ~ 1000, forming a complementary 
cosmological probe. One well-known example of such a test 
is the statistics of strong lensing events. 

This paper investigates the use of strong lensing statis- 
tics as probes of cosmology, and in particular discusses the 
joint constraints that may be obtained on the equation of 
state of the dark energy, P — wp (where in this paper w is 
assumed constant), and on the fractional matter density !5m 
and power spectrum normalisation, parameterised by as, the 
rms linear density fluctuation in spheres of 8h -1 Mpc. Mul- 
tiple images of distant QSOs can be produced by strong 
lensing by intervening mass overdensities such as galaxies 
and clusters of galaxies. The work presented in this paper 
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is focused on wide-separation lensed QSOs, AO <; 5". These 
systems are produced by structures with the mass of clus- 
ters of galaxies, so wide-separation lensed QSOs can probe 
the number and evolution of clusters in a similar way to the 
cluster mass function deduced from X-ray-selected samples 
or statistical measurements of the Sunyaev-Zeldovich effect. 
However, these latter methods depend on emission from, or 
scattering by, the baryons in a cluster, while gravitational 
lensing directly probes its mass. But, despite the existence of 
arcminute-separated multiple images of background gal axies 
JSmail et aljll995t iKneib et alJll9S>d ISahu et alJll99Sl) . the 
detection of wide-separation lensed QSOs has proved to be a 
hard task. Only two systems are known to have separations 
AO > 5": Q0957 + 561, with A9 = 6.26" and redshift z = 
1.41, is the widest-separation confir med lensed QSO known 
JWalsh. Carswell fc Wevmannll 19791) and RX J0921 + 4529, 
with A9 = 6.97" and z = 1.65, remains a wide- 



separation lens candidatej^ur^^^lj 2Q01 I). Other svs- 



tematic lensing searches JOfek et alfl200ltlMaoz et al|l997 



IPhillips et all boOlt IPhillips. Browne fc Wilkinson! 12001 

have failed to find wide-separation lensed QSOs. Ongoing 
redshift surveys such as the Two-Degree Field QSO Red- 
shift Survey (2QZ) and the Sloan Digital Sky Survey (SDSS) 
with 10 4 -10 5 Q SOs might be able t o detect wide separation 
lensed systems. iMiller et alJ i2002f) have reported evidence 
for QSO lensed candidates with arcminute separations se- 
lected from 2QZ, but further observational work is required 
to demonstrate that those candidates are indeed strongly 
lensed systems. 

Stron g lensing statist i cs hav e been calculated for ACDM 
models bvlLi fc Ostrikerl (120021) and for QCDM models by 
ISarbu et alJ (l200lf) . The work presented h e re no t only ex- 
tends the work performed by ISarbu et al.l (l200lT) to wide- 
separation lenses but also uses what we consider to be more 
accurate assumptions about the halo mass function and the 
mass distribution within halos, and investigates more ex- 
plicitly the effect of making differing model assumptions . 
First, we adopt the Sheth & Tormen JSheth et alJl200lf) 
halo mass function, which should more accurately predict 
the mass function of massive clusters. Second, we utilise 
the da rk matter halo concentration prescription of lEke et alJ 
(2001) and take into account the variation of the concentra- 
tion with redshift, mass and cosmology. Here we have explic- 
itly considered the dependence of the concent ration on w, 
Qm an d erg. For comparison wi t h pre vious work. lSarbu et all 
teOOll) used the lBullock et alJ J2001T) relation and only con- 
sidered the mass and redshift dependence of concentration. 
Third, we include th e effects of scatter in halo concentration 
(Bullock et al.ll200 J) wh ich w as neglected bv lLi fc Ostrikerl 
l|2002l) and lSarbu et all J200ll) . These last points are crucial 
as lensing depends very strongly on halo concentration, in 
that more concentrated halos are more efficient lenses. We 
also note that the statistics of strong lensing events depends 
critically on the assumed value of as, and by assuming an 
empirical relationship between as and w the cause of vari- 
ation in lensing statistics with cosmological parameters can 
become obscured. In fact the lSarbu et alJ (12001') results are 
dominated by the assumed as~w relationship. In this paper 
we present results as a function of as, w and Qm so that 
the inter-relationships between these variables may be un- 
derstood. We also show how calculation of the amplification 



bias changes when widely separated multiple images are be- 
ing considered. 

Hence, in section [5] we describe how the lensing prob- 
ability arising from massive clusters may be estimated and 
discuss the elements that play a role in that computation, 
in section |3 we show the results of that analysis and in sec- 
tion 2] we discuss the outlook for being able to constrain 
cosmological parameters from such measurements. 



2 THE LENSING PROBABILITY 

If we define the lens system with a lens at redshift z, a source 
at redshift zs, which is in our case a QSO, the probability 
for the source being lensed with an image separation greater 
than A9 is: 



P(z s , > AO) = 



' (1+zf^dz 
dz 



dti 

— (M,z)dMa le „s(l) 



where ai ens is the lensing cross-section which is specified in 
section 1231 n is the number of lenses per comoving volume, 
Mmin is the minimum mass required to produce an image 
splitting AO, r is the proper cosmological distance to the 
lens and (1 + z) accounts for the fact that n is in units of 
comoving volume. The basic ingredients that take part in the 
calculation of the lensing probability are therefore the halo 
mass function, the background cosmology and the lens cross- 
section which depends critically on the halo profile and on 
the magnification bias. The following subsections will detail 
these ingredients. 



2.1 The Cosmological Model 

We consider spatially flat cosmological models with a dark 
energy equation of state Pq = wpg, where Pq is the pres- 
sure, pq is density and w, the so-called dark energy equa- 
tion of state parameter, is assumed to be constant with 
time. In this case the local energy conservation law (see e.g 
iPeebles fc Ratrall2002T) is pq = —3 (a/a) pq (1 + w), where 
Pq oc a _3 ' 1+ ™- ) and a is the cosmic expansion factor. Dark 
energy has a number of implications in cosmology and some 
of these are of direct relevance for strong lensing. Two 
clear consequences are the effects on the cosmic volume and 
the growth of structure. For a flat model with fixed frac- 
tional matter density Qm (and therefore Qq = 1 — Qm) 
the cosmic volume per unit redshift decreases as w in- 
creases whereas structures start forming earlie r if w increases 
feartelmann et al]l20ol iKlypin et alJl2003l) . For our pur- 
poses we need to compute angular diameter distances. In 
flat dark energy cosmologies the angular diameter distance 
from an object at redshift z\ to an object at redshift 22 is 
given in units of c/Ho by: 



dz 



(1 + z 2 ) J (Qm(1 + zf + Qq(1 + z)3(i+»)) 



Other consequences of dark energy such as its influence on 
the mass power spectrum, the linear and non-linear overden- 
sity for collapse and the dark matter halo concentrations will 
be discussed in later sections. We adopt a value for Hubble's 
constant, Ho = 70kms _1 Mpc -1 . 
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2.2 The Halo Mass Function 

In this paper we consider lensing only by massive galaxy 
clusters, and neglect any contribution from individual galax- 
ies. The results obtained should be valid on image separation 
scales where the calculated lensing probability significantly 
exceeds that due to individual galaxies, and from our analy- 
sis we estimate the minimum scale to be A9 ~ 5". To calcu- 
late the number density and redshift dist ribution of massive 
clusters we adopt the Sheth & Tormen jSheth et aJjl200lh 
mass function. This is a modified version of that obtained 
by the Press- Schechter formalism JPress fc SchechteJll97il 
that fits better the space density of high-mass halos deter- 
mined from N-body simulations. This increase in precision 
is obtained principally by introducing an additional factor 
into the critical linear overdensity for collapse, where that 
factor is adjusted so that the mass function matches the re- 
sults from those simulations. This mass function produces 
results almost indistinguishable from the mass function of 
l.lenkins et al.l (1200 if) . which is also a fit to N-body simu- 
lations. Both fitting formulae predict significantly greater 
numbers of high-mass halos than Press- Schechter. In adopt- 
ing this function, we implicitly have to assume that the func- 
tion is vali d for dark energy mo dels as indicated by the sim- 
ulations of iKlvoin et al.l l|2003ll . The mass function is pri- 
marily a single-parameter modification of Press-Schechter: 
the latter makes no assumption about background cosmol- 
ogy, and already fits remarkably well to results from N-body 
simulations. We can be confident that the mass function pre- 
scription will remain valid for a wide range of cosmological 
models, although we shall argue later that it is important 
to test not only the mass function but also the distribu- 
tion of halo concen tration values by N-body simulation (e.g. 
iKlvpin et aHl2003D . 

We also adopt the linear growth and the line ar matter 
power spectrum for QCDM models obtained by iMa et all 
(1999). The major difference between the power spectra of 
the ACDM and QCDM models is that dark energy may 
cluster spatially on scales corresponding to wavenumbers 
k ;Si 0.01ft Mpc~\ whereas the cosmological constant re- 
mains spatially smooth on all length scales. This signi- 
fies that the QCDM power spectrum on cluster scales has 
the same shape as the ACDM model but has a different 
large-scale shape, and, especially if the normalisation is de- 
termined from large-scal e CMB measur ements, the over- 
all amplitude is chan ged llMa et al.lll999l) . Here we use the 
iBardeen et~afl <1986fl fitting formula for t he ACDM mode l 
transfer function with the modification of ISugivamal (Il995l) 
in order to account for the baryons. The resulting ACDM 
power spectra a r e ext remely close to those computed by 
lEisenstein fc Hul Il999h . 

To avoid varying too many parameters, we fix the pri- 
mordial power spectrum index to be n = 1 and the baryon 
density to be the WMAP best-fit value of Q b = 0.047. In 
fact, varying the primordial index does not change signifi- 
cantly the halo mass function if it is normalised to a partic- 
ular value of as- Modifying the power spectrum to include 
the effects of baryons in fact decreases the lensing probabil- 
ity, somewhat contrary to our initial expectations: this point 
will be discussed again in the context of halo concentrations. 

One critical quantity that varies with w is the lin- 
ear theory critical threshold for collapse of a halo, which 




Figure 1. The Sheth & Tormen mass function and its (lack of) 
sensitivity to the dark energy equation of state parameter w. Thin 
lines are for z = 0.1 and thick lines are for z = 1.0. Solid lines 
correspond to a cosmological constant model with w = —1.0 and 
dashed lines correspond to a quintessence model with w = —0.5. 
The fractional matter density and the power spectrum normali- 
sation were set to H^./ = 0.3 and erg = 0.9 respectively. 



is a key part of the Press-Schechter and Sheth- Tormen 
prescriptions. We adopt the fitti ng function proposed by 
|WeiribCT^^J^arrrionkowski| ll2002li which was modelled af- 
ter |Kitav£mia_&_SjitcJ a 199€a) : in fact the value of the linear 
theory critical threshold for the collapse at the typical red- 
shift of the lens (z ~ 0.5) does not vary considerably, not 
more than approximately 5 per cent, over a wide range of 
Qm and w. Figs and |2] present the Sheth & Tormen mass 
function as a function of halo mass. As is well known, the 
mass function is a steep function of mass, decreases with 
redshift and increases exponentially with as- Fig. |2] shows 
the dependence of the mass function on erg and Qm. Fig. 
plots the mass function at two redshifts, z — 0.1 and z = 1 
for a cosmological model with a cosmological constant and 
for a model with a dark energy equation of state parameter 
w — —0.5. Clearly, the mass function is not very sensitive to 
the dark energy equation of state parameter at z = 0.1 but 
the difference between the two models is greater at higher 
redshifts, and the model with w — —0.5 predicts more halos 
than the cosmologic al constant model, f or a fixed value of 
as, as shown also bv lKlypin et alJ J2003f) . 



2.3 The Lensing Cross Section 

The lensing cross section, which may be defined as the area 
on the lens plane for which multiple imaging occurs, is given 
by: 

aiens =Dl 2 tt /3 A(z, M, (3) d/3 (3) 

Jo 

where Dl is the angular diameter distance from the ob- 
server to the lens, f3 is the angle between the lens and the 
source, j3 cr u{z,M) is the critical angle for multiple imaging 
and A(z, M, f3) is the amplification bias. Note here that the 
amplification bias is part of the integrand and is therefore 
explicitly computed as a function of the angle f3 instead of 
simply calculating the cross section as ai ens — ir(DL(3crit) 2 
as in previous studies. 
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Figure 2. The Sheth & Tormen mass function and its sensitivity 
to Qm and crs. The solid line is for (<jg,Qp„,f) = (0.8,0.3), the 
dashed line for (erg, Qm) = (0.9, 0.3) and the dash-dotted line for 
(<ts,Qm) — (0.9,0.2). Here, the mass function is computed at 
z = and we assume ui = —1.0. 



2.3.1 The Halo Profile 

Following previous lensing probability calculations, we as- 
sume that the lenses have a circularly symmetric surface 
mass density, £(r) = £(|r|), where r is a two-dimensional 
position vector in the lens plane. In this c ase, the lens equa- 
tion reduces to a on e-dimensional form IIMaoz et al.lll997t 
ISchneider et al]ll992l) : 



D s 



(4) 



where /3 is the angle between the lens and the source, Dls 
is the angular diameter distance between the lens and the 
source, Ds is the angular diameter distance between the 
observer and the source, 6 is the angle between the lens and 
the lensed image formed and a is the deflection angle: 

4GM(< r) 



a(r) 



(5) 



where M(< r) is the projected mass enclosed within radius 
r and is given by: 



M(< r) = 2tt / E(r')r'dr' 



(6) 



We discuss the possible effects of departures from smooth, 
circularly-symmetric lenses in section 4. 

The den sity profile of the lenses is mo delled with the 
NFW profile l|Navarro. Frenk fc Whitell997h which appears 
to be a good fit to numerically simulated halos over a wide 
range of masses in various cosmological scenarios. The den- 
sity profile is: 

PS 



p(r) = 



^(1 + ^) 2 



(7) 



where ps and rs are the characteristic density scale and scale 
radius, respectively. This profile is steeper than the singu- 
lar is othermal sphere at large radii but is flatter at smaller 
radii. iMoore et alJ <ll999l) simulations prefer a steeper inner 
profile although the observati onal evidence for the forma- 
tion of radial arcs in clusters JSand et alJl2002l) is favour- 
ing a shallower inner profile. The inner slope is thus an is- 
sue under debate and we will not discuss it here but we 
point out that the lensing probability is rather sensitive to 



its value iLi fc Ostrikerl f2002): lensing probability increases 
as the inner slope steepens because for standard NFW pro- 
files the projected mass density is usually lower than the 
critical density for lensing for all but the innermost radii of 
a cluster. 

The surface mass density of the NFW profile is 
jBartelmannlll99eD : 

2p s r s 



E(x) 



-/(*) 



where x is the radial coordinate in units of scale radius 
i e — . As a result the bend angle takes the form: 



167rp s r 2 g(x) 



a(x) 

v~ x 

where g(x) is ( Bar telmanrj |T996 ) : 



yf(y) 



dy. 



(9) 



(10) 



In order to calculate the lensing cross section it is conve- 
nient to calculate a critical angle for lensing, f3 cr it- Multiple 
imaging occurs if and only if the angle /3 satisfies |/3| ^ Peru 
where /3 C rit is the solution of d(3/d9 = 0. In the case of per- 
fect alignment between the observer, circularly-symmetric 
lens and source, (3 — 0, and \P\ < /3 cr i t there are three im- 
ages formed, for |/3| = f3 C rit two images and for |/3| > f3 cr i t 
only one image is formed. Non-circularly symmetric lenses 
can produce more than three images. A key question for the 
lensing cross section is how to calculate the inner parame- 
ters p a and r a and their variation with mass, redshift and 
cosmological model. These parameters can be related to the 
virial parameters of the halo through the halo concentration 
parameter. The latter is basically a measurement of the rel- 
ative size of the core with respect to the virial radius. Here 
we adopt the approach of calculating the halo concentration 
from prescriptions based on N-body simulations. 



2.3.2 Dark Matter Halo Concentrations 

Formally, the mass of the NFW profile diverges, so to define 
the mass we adopt the standard formalism of calculating the 
mass within the virial radius Rvi r of a sphere with density 
A V i r times the critical density p cr it'- 



M vir = 4tt 



pr 2 dr = A-Kp s r 3 s f{c vi 



47T 

^ ^ vir per i t 



(ii) 



and we define the halo 



where f(c vir ) = ln(l + c vir ) — j 1 
concentration parameter as c V i r — R ^* r . A„i r is the non- 
linear overdensity at virialisation. With these definitions we 
can compute the scale radius r s = R V i r /c V i r and the char- 
acteristic density ps- 

„3 



As 



pcrit{z) 



(12) 



In the work that follows we take into account the variation 
of the mean halo concentration with redshift, halo mass and 
cosmology, and in particular with w. The assumptions on 
these dependencies turn out to m ake a significan t differ ence 
to the overall lensing probabilities: ILi fc Ostrikerl (l2002f) had 
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Figure 3. Dark matter halo concentrations function for four dif- 
ferent combinations of redshift and dark energy equation of state 
parameter (z,w). The fractional matter density and the power 
spectrum normalisation were set to Qm =0.3 and erg = 0.9 re- 
spectively. Thin lines are for z = 0.1, thick lines for z = 1.0. Solid 
lines correspond to a ACDM model (ui = —1.0) and dashed lines 
correspond to a quintessence model with w = —0.5. 

no depe ndence of concentr ation on mass in their ACDM 
analysis. ISarbu et alJ i200 ll) included the redshift and mass 
variation in their QCDM studies but did not include any 
dependence on w, assuming instead the values appropri- 
ate for a ACDM cosmology. To evaluate those dependen- 
cies w e adopt the concentration prescription of lEke et alJ 
teOOll) . That prescription is based around the premise that 
halo concentration increases as the redshift of formation of 
a dark halo increases, a reflection of the higher matter den- 
sity at earlier epochs. A characteristic redshift of formation 
that varies with the shape of the power spectrum is cal- 
culated and a value for a single parameter can be found 
which yields concentration values that match well the re- 
sults from N-body simulations for a wide range of cosmolog- 
ical models, and has been tested against the original NFW 
results for masses up to 10 15 Mq. Concentrations deduced 
in this way seem to be generally in broad agreement with 
values measured from fitting to the X-ray profi les of clus- 
ters o ver the range 10 14 M vir ^ 1O 16 M (e.g. lWu fc Xud 
2002). It does appear that the dependence on mass in those 
samples may be steeper than predicted, but X-ray selec- 
tion of clusters is likely to have some effect on the observed 
relationship, and at high masses, which are most impor- 
tant for wide-separation lensing, comparable concentration 
values are obtai ned from the measurements and from the 
lEke et afl (1200 J) prescription. 

In this paper we extend the lEke et alJ (1200 ll) prescrip- 
tion in a natural way to accommodate dark energy mod- 
els working under the assumption that the algorithms de- 
rived from numerical simulations in ACDM models are also 
valid in QC DM models. We have used the fitting func- 
tion of IWeinberg fc Ka mionkowskil i2002l) . modelled after 
lKitavam^^^uton^96^ for the non-linear overdensity at 
virialisation, A„j r , which increases with w. This arises be- 
cause structures start forming earlier in models with a 
higher w, the mean energy in a collapsing object is larger 
and hence a higher overdensity fo r virialisation is required 
iWeinberg fc Kamionkowskil |20O2) . Figures [3] and 2] show 
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Figure 4. Dark matter halo concentrations for three (cts,Qm) 
combinations. The solid line is for (<Tg,f2^f) = (0.8,0.3), the 
dashed line for (ffg,!!^) = (0.9,0.3) and the dash-dotted line for 
(<rg, £Im) = (0.9,0.2). Concentrations were computed at z = 0. 
and assumed a ACDM model. 

the concentrations obtained for differing cosmological mod- 
els. The concentrations decrease with the mass of the halo 
since high halo masses are formed later in structure forma- 
tion models. In quintessence models structures form earlier 
than in ACDM, resulting in higher concentrations for a given 
ma ss. This trend is evident in the plot of Fig. |2l Although 
the lEke et al.l i200ll) prescription has not yet been exhaus- 
tively tested for w models, comparison with the dark en- 
ergy simulations of iKlypin et alJ (120031) show good general 
agreement, with higher w models exhibiting higher concen- 
trations for a given mass. Fig. 2] presents concentrations for 
three combinations of (cts,Qm)- An increase in <rg makes the 
halo concentrations higher and a decrease in !5m lowers the 
concentrations. 

One key additional factor to consider is the scatter in 
halo concentration values about the mean relations com- 
puted above. The numerical simulations show that there is a 
scatter in the concentration values which may be associated 
with a spread in col l apse epochs and /or f ormation histories 
jBullock et al.ll200ltlWechsler et alJl2002T) . This scatter has 
a dramatic effect on the wide-separation lensing probability 
because clusters with NFW profiles and moderate concen- 
tration values have surface mass densities which only rise 
above the critical density for lensing at small radii. Even a 
30 percent increase in concentration can bring the mass den- 
sity above the critical value over a substantially larger range 
of radii, thereby greatly increasing the lensing probability 
at large separations. So giving a minority of the cluster halo 
population larger concentrations than the mean produces 
a significant increase in lensing probability, as discussed in 
section 3. The scatter also reduces the differences between 
the various models however, as also discussed later in this 
paper. The amount of scatter measured from N-body simula- 
tions depends critically on which halos are chosen for study: 
relaxed halos are better fit by NFW profiles and have a 
moderate scatter. iBullock et alJ (1200 ll) argue that the scat- 
ter appears to have log- normal distribution with a c = 0.18 
and we incorporate this distribution into the com putation of 
the le nsing probability. Finally we note that the lEke et alJ 
teOOll) prescription has been tested in a range of power spec- 
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tra shapes. Here we find that introducing Qt = 0.047 lowers 
the mean halo concentrations and therefore results in a de- 
crease of the lensing probability. 



2.3.3 The Amplification Bias 

The effect of the amplification bias on the lensing cross sec- 
tion, A(z,M,f3), is to boost the lensing probability because 
sources that are intrinsically dimmer than the flux limit of 
the survey fu m are brought into the sample by the mag- 
nification. If the source QSOs have a single power-law flux 
distribution with redshift-independent slope a, N zs (> /) = 
f~ a then the amplification bias is given by: 



A(z,M,(3) 



N zs (> f Hm /fi(z,M,l3)) 

Nz s (> flim) 

= v a {z,M,(3) 



(13) 



where fi is the magnification factor which for circularly sym- 
metric lenses is given by: 

dO 



( D LS da(6) 
V Ds d9 



- f 



(14) 



In order to calculate the magnification factor we need to 
solve the lens equation and obtain the three (or two) solu- 
tions. In surveys for wide-separation lensed systems, where 
the image separation is larger than the survey resolution, 
the multiple images are separately detected and we need 
to compute the probability that at least two images will 
appear brighter than the survey flux selection limit. Hence 
the relevant quantity to compute is the amplification of the 
second-brightest of the multiple images. This approach dif- 
fers from that commonly taken in which multiple images are 
assumed to be unresolved and the amplification bias is cal- 
culated from the sum of the fluxes of the lensed images. In 
addition, in this paper we compute both the magnification 
factor and the minimum mass required to produce a given 
image splitting as a function of the angle (3, in contrast to 
the usually-followed procedure of calculating the minimum 
mass for lensing assuming perfect alignment between the 
observer, the lens and the source. 



3 RESULTS 

Having explored the sensitivity of the mass function and the 
dark matter halo concentrations to the dark energy equation 
of state parameter as well as to the fractional matter density 
and the power spectrum normalisation, we now investigate 
the influence of the same cosmological parameters on the re- 
maining lensing ingredients, such as the lensing cross section 
and the minimum mass for lensing. 

We first look at the lensing cross section: in order to 
separate the various effects we do not at this stage ap- 
ply the magnification bias, and hence effectively calculate 
o~iens = Tr(DLf3crit) 2 ■ Fig. [SJshows that the lensing cross sec- 
tion increases monotonically in quintessence universes for 
fixed ag, and the amplitude of the cross section is larger 
if the power spectrum normalisation is raised. This arises 
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Figure 5. The lensing cross section as a function of w for a halo 
of mass Mualo = 10 14 Af Q and redshift z = 0.3. The QSO is 
at redshift zqso = 1-5- Curves are for three different (ag,Q.M) 
combinations. Solid curves are for Qm = 0.3 and erg = 0.8, dashed 
for Ojvf = 0.3 and ag = 0.9 and dash-dotted curves for Qm = 0.2 
and ag = 0.9. No amplification bias has been applied. Thin lines 
show the lensing cross section cr; ens = 7r(D^/3 cr i t ) 2 . Thick lines 
include the additional effect of introducing the scatter in the halo 
concentrations. 
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Figure 6. The minimum mass required to lens a QSO at redshift 
z QSO = 1-5 with image separation A8 = fO". The lens is at 



redshift z^ e 



0.3 and it is assumed that there is close to 



perfect alignment, /3 ss 0. The minimum mass is plotted as a 
function of w for three pairs of (ag, Qm)- The thin solid line is for 
(ag,Q M ) = (0.8,0.3), the thick solid is for (ag,Q M ) = (0.9,0.3) 
and the dashed line is for (ag,QM) = (0.9,0.2). 



because the concentrations are larger in those cases which 
in turn make lensing more efficient. If the matter content 
is decreased from Qm — 0.3 to Q m = 0.2 the cross sec- 
tion is lowered, although the difference between the models 
decreases for higher values of w. This occurs because halo 
concentrations are in fact lower for lower values of !Jm (see 
Fig. Also shown is the additional effect of including a 
scatter in the value of the concentrations, in particular by 
assuming the log-normal distribution with a c = 0.18. In- 
cluding the scatter significantly increases the lensing cross 
section: for example, in a model with (o^Om) = (0.9,0.3) 
and w = — 1 its value is raised by a factor of approximately 
one hundred, as the tail of high concentration halos domi- 
nates the statistics. We should expect that for large scatter 
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Figure 7. The lensing probability as a function of image separa- 
tion A9 for four different values of the slope of the QSO number- 
counts. The solid line shows a = 0, corresponding to no amplifica- 
tion bias. The fractional matter density and the power spectrum 
normalisation were set to Qm = 0.3 and rrg = 0.9 respectively. 
The dark energy equation of state parameter was set to»= —1. 



in halo concentration the differences between cosmological 
models should be reduced (assuming the scatter itself is in- 
dependent of cosmology) and this is reflected also in Fig. |S] 
Concentration scat ter was not included in pr evious work by 
ISarbu et all <200ll) and lLi fc Ostrikerl J2002I) . 

The other key ingredient on the lensing calculation is 
the minimum mass required to produce multiple imaging 
and this is included in formula as the lower limit of the 
mass integral. This is generally calculated assuming perfect 
alignment between the lens, observer and the source, that 
is to say with j3 = 0. For the lensing probabilities here we 
actually compute the minimum mass as a function of the 
angle (3 but in order to analyse how its value varies with 
cosmological parameters we assume (3 = 0. Fig. HO shows 
that the minimum mass required to lens a QSO at redshift 
zqso = 1-5 with an image separation A6 = 10" and a lens 
redshift zi ens = 0.3 clearly decreases with w and a%. For 
(7s = 0.9 the minimum mass required is greater for a universe 
with a lower fractional matter density. 

We can now look at the probability distribution of im- 
age separations fFieslTland llOII . There is a clear and steep 
fall of the probability toward higher image separations which 
is a consequence of a mass function that falls off quite rapidly 
combined with a higher required minimum mass for lens- 
ing events with larger separations. Here we underline the 
sensitivity of the lensing probability to the slope of the 
source number-counts. Fig.|7|shows the expected behaviour: 
a steeper slope favours higher lensing rates. Note that if mag- 
nification bias is not accounted for then the lensing proba- 
bility is reduced by a factor of about ten. At faint radio 
and optical selection levels, with B gi 20, the slope of the 
QSO number-counts are significantly flatter than Euclidean, 
although because we need to know the slope at flux levels 
fainter than the unlensed limit for the survey being searched 
for lensed QSOs the actual value of a is often uncertain. The 
observed apparent slope for the 2QZ optical QSO survey 
has value a m 0.8 at faint magnitudes, although this value 
does not take account of sample incompleteness. Hereafter 
we shall assume a slope a — 1: to first order the change in 




^QSO 

Figure 8. The lensing probability as a function of source redshift, 
zqso, for three different image separations: A8 = 5" (top line), 
A8 = 10" (middle line) and A8 = 30" (bottom line), assuming 
Q.M = 0.3, erg = 0.9 and w = — 1. 
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Figure 9. The lensing probability as a function of concentra- 
tion scatter, <r c , for three different image separations: A8 = 5" 
(top line), A8 = 10" (middle line) and A8 = 30" (bottom line), 



assuming !1m = 0.3, as = 0.9 and w - 



-1. 



lensing probability resulting from assuming a different value 
for a may be estimated from Fig. Q 

Assumptions about the typical source redshift and the 
distribution of source redshifts also affect the lensing proba- 
bility. Fig. |H] shows that the lensing optical depth is a steep 
function of source redshift: if a source is moved from red- 
shift z — 1.5 to z = 2.5, the lensing probability is increased 
by approximately a factor of three. This behaviour is also 
found in estimates of the lensing optic al depth for giant arcs 
JWambseanss. Bode fc OstrikerbooH) . Fig. 8 shows that ap- 
plying an accurate estimate of the median source redshift is 
essential. But we also expect that a distribution of source 
redshifts about the median would increase the lensing prob- 
ability, and indeed if w e adopt the redshi ft distribution of 
the 2QZ 10K catalogue llCroom et all200ll) we find that the 
lensing probability is increased by approximately 10%. 

Uncertainties in halo concentrations which can be ex- 
pressed by <7 C can also alter the lensing probability. This 
is shown in Fig. [5] if the scatter in halo concentrations is 
changed from 0.18 (the value assumed in our analysis) to 
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Figure 10. The lensing probability as a function of image sep- 
aration A8 for four different values of id in a flat universe with 
<T8 = 0.9 and Qm = 0.3. 



0.4 the lensing probability is increased by a factor of six. 
Therefore, an accurate value of the scatter in halo concen- 
trations is fundamental. 

The effect of increasing w > — 1 on the probability dis- 
tribution of image separations is to boost the integrated lens- 
ing probability. The difference between the models increases 
for larger w (Fig. llOt . This is partly owing to halo concentra- 
tions being larger in quintessence models, making the lensing 
cross section higher, and partly owing to the minimum mass 
required for multiple imaging decreasing with w (Fig.|SJ: the 
mass function of galaxy clusters falls sufficiently steeply that 
the small decrease in mass has a significant effect. 

Fig. llll explores the sensitivity of the lensing probability 
to Qm and ag for differing minimum values of separation, 
A8 = 10" and AO = 30". The latter probability is about 
ten times smaller than the former. We confirm the expected 
increase of probability with w as in Fig. 1101 and note that 
the probability also increases with ag. Indeed the lensing 
cross section increases with <rg and the minimum mass for 
lensing also decreases with this cosmological parameter. In 
addition, there is an exponential rise of the mass function 
with its value. The overall effect is therefore an increase of 
the lensing probability. When we lower Qm from 0.3 to 0.2 
the probability decreases. 

There are evident cosmological parameter degeneracies 
and those are clearly shown on the contour plots provided in 
Figs 1121 and 1131 One has to bring in other independent cos- 
mological constraints in order to break the degeneracies. In 
Fig. ll2l we overplot the best-fit r esult from the abundance o f 
rich galaxy clust e rs obt ained bv lWang fc Steinhardtl jl99ct) 
and lLokas et alJ i2003t) as well as the constraint obtained 
from th e COBE measurement of the cosmic m icrowave back- 
ground iMa et alJll999l:lBunn fc Whitdll997l) . In Fig. EH we 
show the constraints from a st udy of the present -day num- 
ber density of galaxy clusters jAllen et al.|| 2003|) and from 
a cosmic shear measurem ent jHoeltstoa J _Yee anriGladdersI 
l2002l) . |Jarvis et alJ i2002T) and lPierpaoli et alJ feOOll) review 
the variation found in these measurements from studies per- 
formed by other groups. T he WMAP best-fit result is also 
shown dSpergel et al.ll2005t . 

The ability of surveys to actually discriminate between 
cosmological models will depend on survey size and on the 
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Figure 11. The probability of a QSO at redshift zqso = 1-5 
being lensed with an image separation bigger than A8 in a 
quintessence model. The sensitivity of the probability to erg and 
Qm is also shown. Solid lines are for Qm = 0.3 and erg = 0.8. 
Dashed lines are for Qm = 0-3 an d CT 8 = 0.9 and dash-dotted 
lines are for Qm = 0.2 and erg = 0.9. Thick lines correspond to a 
minimum image separation of A8 = 30" and thin lines correspond 
to Ad = 10". 
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Figure 12. Contours of probability P(zq SO = 1.5, AS > 10") 
in the w — erg plane for a universe with fractional matter density 
Qm = 0.3. The dashed and dotte d lines show the constraints from 
cluster abundance s as obtained by Wang & Stcinhardt ( 1998]) and 
lLokas et alJ 120031) . respectively. The dashed-dotted l ine shows 
the r esult of COBE power spectrum normalisation jMa et alJ 
1999). 



numbers of detected lenses of course. As an illustration, we 
assume a survey of 10 5 sources containing four lensed sources 
with A8 > 5", which corresponds to the model prediction 
for a flat cosmological model with w — —1, erg = 0.9 and 
Qm = 0.3, and the source redshif t distribution of the 2QZ 
10K catalogue JCroom et alJl200ll) . The likelihood function 
chosen depends on the number of lenses as 



L 



A e 



Nl 



where A is the expected number of lensed systems and 
is the actual number observed; we do not here consider in- 
formation contained in the distribution of angular separa- 
tions. Fig. 1141 shows the 68.3 and 95.4 per cent confidence 
level contours in the erg — Q m plane. Fig. H5l shows the con- 
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Figure 13. Contours of probability P(zqso = 1-5, A8 > 10") in 
the as — Ojvj plane for a universe dominated by a cosmological 
constant (w = —1). The area enclosed by the dashed lines rep- 
resents the constraints from the study of the local cluster abun- 
dance I Allen et af]|2003l) . The area enclosed by the dot-dashed 
lines shows the constraints obtained from the c osmic shear mea- 
surement of lHoekstra. Ye e and Gladdersj l| 2002fl . The dotted line 
is the WMAP relation iSpergel et all2003l) . 

straints on w when one assumes Qm = 0.3 and as = 0.9 
(marginalising over the current uncertainty in as would sig- 
nificantly broaden the likelihood function). The confidence 
intervals are somewhat broader than can be attained from 
current generations of cosmic microwave background, weak 
lensing measurements or deductions from cluster abundance 
measurements, but it is important to note that these meth- 
ods are all complementary tests of the cosmological model. 
The CMB measurements probe the Universe at z ~ 1000; 
cluster abundance determinations are dependent on correct 
modelling of thermal emission from baryons in cluster po- 
tential wells; cosmic shear measures the mean fluctuations 
in matter density on Mpc scales; the abundance of strong 
lensing events measures the extreme high-mass tail of the 
mass function, with a small number of lensed systems being 
generated by the most extreme objects. In the final analy- 
sis the abundance of strong lensing events may tell us more 
about the formation of the most massive structures in the 
universe than about the values of cosmological parameters. 



4 DISCUSSION AND CONCLUSIONS 

We have calculated cluster lensing probabilities of QSOs in 
dark energy models. The lensing probabilities have been ex- 
plored over a wide range of angular separations and we have 
taken into account the cosmology, mass and redshift depen- 
dence of halo concentration, while assuming that the algo- 
rithms derived from dark matter halo simulations in ACDM 
models remain valid in QCDM models. Ultimately, these 
algorithms should be tested against numerical simulations 
in dark energy cosmologies, but simulations to date indi- 
cate that these proscriptions are at least approximately cor- 
rect IIKlvpin et al. 2003). What issues remain that could still 
have an important effect on the predictions of this type of 
modelling? At this stage we may be reasonably confident 
that the Sheth & Tormen mass function is sufficiently well 
tested that it can be considered a robust prescription, but 



Figure 14. The 68.3 (solid line) and 95.4 (dashed line) per cent 
confidence level contours on as and Qm (assuming w = — 1.) 
obtained for a survey of 10 5 sources, four of which are lensed 
with A6» > 5". 
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Figure 15. The 68.3 (w < -0.66) and 95.4 (w < -0.51) per cent 
confidence level likelihood constraints on w (assuming Qm = 0.3 
and as = 0.9) obtained for a survey of 10 5 sources, four of which 
are lensed with A8 > 5". 

we should be aware that modifications to the fluctuations 
giving rise to cosmological structure would lead to differing 
mass functions. As one example, non-Gaussian initial per- 
turbations significantly increase the mass function at high 
redshifts (iMatarrese et al.ll200ol ; iMathis et afll2003l) and we 
should expect halos of a given mass to have formed at earlier 
epochs, and hence have higher values of halo concentration. 

There is currently some uncertainty about the univer- 
sality of the NFW profile and about its applicability to lens- 
ing studies. Merely by introducing the known scatter in halo 
concentration we see a significantly non-linear effect on lens- 
ing probability. If, in turn, this scatter is greater then the 
value adopted in this paper (a c = 0.18) then that will have 
a great effect in the lensing probability. As we have seen in 
Fig. |5| if the scatter in halo concentrations is changed from 
0.18 (the value assumed in our analysis) to 0.4, the lensing 
probability is increased by a factor of six. Variations in the 
inner slope of the halo profile would also have a significant 
effect iLi fc Ostrikedl2002l : lQgurill2003al) . There have been 
claims based on N-body simulations that th e inner slope 
may be steeper than the standard NFW valu e llM oore et all 
1999), or that its value steepens with mass (lRicottill200*2T) 7 
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and there is evidence that some clusters at least m ay have 
slopes that are flatter than NFW JSand et alJl2002l) . Even if 
there is not a systematic departure from the canonical NFW 
profile, scatter about that relation could lead to increase in 
lensing probability in a similar way to that introduced by 
the scatter in halo concentration. Other effects to be taken 
into account include the effect of substructure, especially in 
recently-formed, non-relaxed clusters: this has already been 
show n to be important for t he statistics of extreme giant 
arcs jMeneghetti et al . 2003a). Finally, the effect of individ- 
ual galaxies in a cluster may be noticeable on the smallest 
angular scales discussed in this paper, although present in- 
dications are that this is a relative ly small effect in giant arc 
statistics jMeneehetti et al-l feoOO) and, likewise, the inclu- 
sion of a massive cD galaxy in the centre of the cluster does 
not seem to si gnificantly enhance the lensing cross section 
for giant arcs iMeneehetti et alJl20 03b'l. 

We have also adopted a circularly symmetric lens pro- 
file, whereas even relaxed individual clusters are likely to 
be triaxial. Again, this effect has been shown to be po- 
tentially significant for the statistics of extreme giant arcs 
(Bartelmann et al. 2002; Oeuri 2003b), although it is not yet 
clear whether there is any significant effect on the statistics 
of less extreme lensing events. 

In applying these models to actual surveys, we should 
of course also take care to ensure that we are using accu- 
rate estimates for the slope of the QSO number-counts and 
that we take account of the redshift distribution of QSOs. In 
particular, an accurate estimation of the median redshift of 
the sources is essential as the lensing probability is a steep 
function of source redshift. For optical QSO surveys this is 
known exactly and therefore does not constitute a source of 
uncertainty. However, for radio surveys the median redshift 
of the sources can be quite uncertain and therefore provide 
a systematic error in the lensing probability: if the median 
redshift is increased from z — 1 to z = 1.5 the lensing prob- 
ability is increased by a factor about three (Fig. [SJ. 

With these caveats, we can use the models produced 
here to understand what constraints, if any, may be placed 
on cosmological models as a result of either detecting or not 
detecting wide-separation lensed QSOs in surveys such as 
SDSS and 2QZ. Within the cosmic concordance model, a 
flat model with Q M = 0.3 and a WMAP value of ag = 0.9 
we find that the lensing probability increases with the dark 
energy equation of state parameter and attribute this effect 
to the fact that the concentrations of dark matter halos are 
higher in quintessence models and therefore lenses become 
more efficient. This effect combined with a mass function 
that is very steep and a lower minimum mass required for 
multiple imaging makes more halos available for strong lens- 
ing. We also note that the ag-flM degeneracy is equivalent 
to the degeneracy found in other methods such as cluster 
abundances, cosmic shear measurements and the WMAP 
analysis (Figs. Upl and 1140 , 

In considering whether lensing statistics can distinguish 
between values of w, we see that at present the uncertainty 
in the calibration of erg dominates the conclusions reached. 
Cluster-normalised models predict that lensing probability 
is not a sensitive indicator of the value of w, whereas COBE- 
normalised values of ag indicate that as w increases, w > — 1, 
the expected number of lensing events decreases. Further 
work is needed to tie down the ag normalisation in w mod- 



els. Of course, we should not expect cluster- and CMB- 
determined normalisations to agree except at a single point 
in parameter space that corresponds to the values for the 
universe we actually inhabit, but at present we cannot be 
confident that such a point can be identified. Nonetheless, 
the determination of the absolute number of lenses at red- 
shifts ziz3 should provide a good test of the ACDM model, 
independently of, and in direct contrast with, cosmic mi- 
crowave background anisotropy measurements at z ~ 1000. 
Alter native methods, such as 3D w eak lensing ([H eavensl 
2003) and weak lensing tomography jjain fe Tavlorll2'o03TK 
although at present with current surveys are not able to 
constrain w, potentially offer a more robust probe with an 
estimated accuracy for w better th an 5% for surveys su ch 
as the one proposed with the LSST (Tain fc Tavlorll2003ft . 

For a flat cosmological constant model with Qm = 0.3, 
a WMAP value of ag = 0.9 and t he source redshift d istri- 
bution of the 2QZ 10K catalogue ijCroom et all 1200 lft . the 
modelling predicts probabilities of 4.0 x 10~ 5 , 2.5 x 10~ 5 
and 4.7 x 10 -6 for lensing events with separations larger 
than 5", 10" and 30" respectively. This indicates that a sur- 
vey with 10 5 sources will find approximately 4 lensed QSOs 
with separations greater than 5", 2 lensed QSOs with sepa- 
rations greater than 10" and none with separations greater 
than 30". For surveys such as the 2QZ and SDSS the survey 
selection effects will somewhat alter these numbers. If no 
wide separation lensed QSOs are detected that would rule 
out models with high ag. Conversely, if large numbers are 
found, or if any lensed systems significantly larger than 1' 
are confirmed, we would need a significant modification to 
the model presented here, including either significant effects 
from the presence of substructure or perhaps requiring as- 
sumptions in the standard model, such as Gaussianity, to 
be relaxed. At this stage it would be premature to specu- 
late further. Challenges in the immediate future, then, are 
to test the expected levels of lensing probabilities against 
numerical simulations, both to take account of substructure 
and cluster shape, and to test the predictions in the cases 
of dark energy models. Joint comparison between multiple 
image probabilities and the statistics of extreme giant arcs 
would be a good simultaneous test. 
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